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Abstract 

We examine the magnetic response of free clusters considering the spin di- 
rection and the cluster orientation as the only active degrees of freedom. 
The average magnetization in small fields approaches the Langevin value for 
paramagnets, depending on the degree to which the Hamiltonian preserves 
symmetries. Superparamagnetic behavior is not achievable within models 
considering only these degrees of freedom. 
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I. INTRODUCTION 



Experimental studies of magnetism of free atomic clusters have shown a number of 
interesting phenomena Susceptibilities measured in Stern-Gerlach experiments 

can be large, with an irregular dependence on the number of atoms in the cluster. High 
susceptibilities are found not only in normally ferromagnetic elements but also in an element, 
rhodium |Q , which is not ferromagnetic in the bulk. The character of the susceptibility is 
also variable. In many cases it is described phenomenologically as superparamagnetism, 
which means that besides a large susceptibility the clusters behave uniformly, independent 
of initial state within the thermal ensemble. The clusters may also respond with a large 
variance in their magnetization, with some clusters even having magnetic moments opposite 
the field. Of course, the limiting behavior of this kind is the magnetization distribution of 
an isolated atom, proportional to the spin projection along the magnetic field. 

The general theoretical problem we address is whether the observed behavior of free 
magnetic clusters can be understood in the framework of Hamiltonians that treat only rota- 
tional and spin degrees of freedom, neglecting vibrations. Vibrational coupling is of course 
important for the response of warm clusters [Q, but the question of whether cold clusters can 
show superparamagnetism is still open. This paper is an extension of earlier work which con- 
sidered a very simplified Hamiltonian with only rotational degrees of freedom active |^,|Tn . 
Large average susceptibilities were found to be possible, but not superparamagnetism. 

We shall calculate the cluster magnetism in a quantum mechanical representation, al- 
though the angular momenta are high enough that classical methods should suffice. The 
rotor is represented in the usual basis {J KM), with M the component of angular momentum 
about the z axis and K the component about the intrinsic third axis. The rotor is coupled 
to spin states in the basis {SSz)- Taking the field direction along the z axis, the component 
of angular momentum about that axis, M + Sz-, is conserved. In the earlier studies, the 
simplified Hamiltonian conserved K as well. 

The theoretical magnetic response associated with a particular Hamiltonian is most 
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conveniently discussed by comparing with Langevin model of paramagnetism. Here the spin 
is coupled to a nonspecific heat bath to give the following formula for the magnetization in 
an external field, 

, ^ _ I^pdn, fi,exp{fi,B/kT) 



(coth(x)--) (2) 



Here /io is the intrinsic magnetic moment and the parameter x is the ratio of magnetic 
energy to thermal energy, x = fi^B / kT. This formula is sometimes used to deduce intrinsic 
magnetic moments from the observed defiections in the Stern-Gerlach apparatus. At low 
field strength the Langevin formula has the limiting behavior, 

M = ^. (3) 

Ho 3 

The condition x -C 1 is in fact satisfied for the conditions of the cited experiments, which 
involved clusters of a few tens or hundreds of atoms. 

Given the Hamiltonian, the magnetic response is calculated as follows. First, the Hamil- 
tonian is diagonalized for a range of magnetic field strengths to produce a spectrum of 
eigenenergies Ei{B). The magnetic moments of the states are then calculated using the 
formula 

^l,{B) = m^^zVB) = (4) 

Since eq. (1) does not require the eigenf unctions, a substantial savings of computation time 
is possible with respect to a full diagonahzation. The distribution of eigenvalues is given by 
the expression 

Here is the dimensionality of the Hamiltonian and Pi are the probabilities of the states. 
We take the adiabatic limit by assigning the probabilities as the Boltzmann factors of the 
levels in the absence of a magnetic field. 



Pi =exp(-^,(0)/fcT). 



(6) 



In applying this, the levels are always ordered and indexed by energy. Thus the lowest level 
in zero field sets the Boltzmann factor for the lowest level at other field strengths. This 
assumes that all level crossings are avoided crossings. It does appear for the Hamiltonians 
we study that the only true level crossings are those imposed by symmetry (M and possible 
K). We do not address the question of whether the field changes take place slowly enough 
for the adiabatic assumption to be valid at the avoided level crossings. 



We consider a Hamiltonian model in which the spin orientation is coupled to an intrinsic 
axis of the cluster and to an applied field B. The Hamiltonian can be written as 



where Jj's are the moments of inertia. The second term is the anisotropy energy, the coupling 
of the spin to the intrinsic axis of the cluster. We assume a uniaxial anisotropy having the 
form of an L = 2 tensor, which would apply to clusters of low symmetry. A Hamiltonain 
of this form was also considered in ref. [0. The usual anisotropy energy applying to cubic 
crystals has the form of an L = 4 tensor. 

The dimensionalities of the state vector spaces for cluster spin-rotor wave functions 
require us to make additional simplifications. The angular momenta present in these clusters 
are of the order of 1000 for a cluster of 100 iron atoms at room temperature. It is impractical 
from a numerical standpoint to deal with ensembles having such large angular momenta. 
Since the quantization of angular momentum should not play a role at large values, we 
believe it is acceptable to use a smaller basis, scaling the parameters of the Hamiltonian and 
the temperature appropriately. Even with a smaller basis, the full Hamiltonian still gives 
large matrices, so we consider as well a limiting case where the spin is not an independent 
dynamical variable, but is fixed to the cluster axis, the "locked-spin" model. 



II. HAMILTONIAN 
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We further assume in the locked-spin model that the cluster has two equal moments of 
inertia, which we take as the 1- and 2- axes, Ii — h- Then in the usual \ JKM) basis K is 
conserved by the rotational Hamiltonian. Its nonvanishing matrix elements are the diagonal 
elements, given by 

A fully asymmetric rotor requires an additional term in the Hamiltonian with matrix ele- 
ments proportional to 



{JK + 2M\Jt - J^\JKM) = ^J{J{J + 1) - K{K + 1))(J(J + 1) - {K + 1){K + 2)) (9) 

We will briefly consider this term later. 

The rotor is coupled to the magnetic field through the spin, whose orientation with 
respect to the intrinsic 3-axis is denoted by 9. The matrix elements of the magnetic field 
term in the Hamiltonian are diagonal in M but not necessarily in J or K. They are given 

by 



{JKM\i^,B\J'K'M') = 6m,m'I^oB^'^j^ {J'MW\JM) 

xl^cose{J'K'lO\JK) + ^[{J'K'l-l\JK)-{J'K'll\JK)]j . (10) 

The more general Hamiltonian allowing intermediate coupling of the spin is constructed 
in the \ JKMSSz) basis. The magnetic energy is trivial in this basis. It depends only on Szi 
with the matrix element {Sz\^qB\Sz) — /IqBSz. The anisotropy energy is given by 



{JKMSSz\j:mYLs2m\J'K'M'SS'z) = (-l)25+^^-^.yX / (^^+l)(^+;)(^,^+3)5(25-l) 



Sk,k'Sm+s.,m'+s'M M 2 (M'-M)l J' M'){J K 2 0| J' K){S 2 {M-M')\S S^) 

(11) 
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III. RESULTS: LOCKED SPIN MODEL 



The locked-spin model for a spherical rotor was treated in ref. [^,20,10,12]. In this model, 
K is a conserved quantum number if the intrinsic 3-axis is taken along the spin direction. The 
previous work showed that the adiabatic ensemble produces an average susceptibility in low 
fields of 2/3 of the Langevin value. The distribution P{n) is very broad, however, far from 
the idealized superparamagnet. A weak point of the model is its assumed K conservation: 
one could argue K conservation is too restrictive to allow adequate mixing of the spin with 
the other degrees of freedom. We shall find that breaking K does in fact have a significant 
effect on the average magnetization, but not on the spread. 

We first remark that changing the Hamiltonian by making the rotor deformed has neg- 
ligible effect on the distribution, in the absence of K mixing. When K is conserved, the 
deformed rotor has the same wave function as the spherical rotor. The magnetic moments of 
the states are thus the same, only the occupation probabilities differ. The thermal ensemble 
would then favor states a somewhat different distribution of K values, but the overall effect 
is small. This is illustrated with the moment distributions for the deformed symmetric rotor 
shown in Fig. 1. The curves were calculated with an ensemble having a cutoff angular 
momentum of 40 for the rotor. The temperature was fixed at kT = 100 in units of If ^ and 
the magnetic field strength has the same value, fi^B = 100/f ^. To display the distribution, 
eq. d^) was smeared with a Gaussian of the form exp(— (/x — /ij)^/(A/i)^) with A/i = 0.05/io. 
The different curves show results for Hamiltonians with moment of inertia ratios ranging 
from 1/2 to 2. One sees that there is hardly any change in the distribution of magnetic 
moments. The average magnetization changes by less than 10% over this extreme range of 
deformations. 

Without discussing in detail the i^-mixing in the Hamiltonian, it is possible to see in a 
crude model how it affects the magnetization. We compute the energy with a Hamiltonian 
that conserves K, but combine the ensembles for different K (but the same M) into a 
common ensemble in eq. (Bl). This assumes that the K-mixing makes all level crossings 
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avoided crossings, but that it is not strong enough to affect the states between crossing 
points. This ensemble produces a higher average magnetization, as might be expected: the 
crossing of states of different K allows the probability to jump to states of larger magnetic 
moment. In Fig. 2 we show the average magnetization of this mixed- ensemble compared 
with the Langevin and with the fixed-i^ locked-spin model. One sees that the mixed-ii' 
ensemble has a susceptibility close to the Langevin model for field strengths below x = 0.2. 
However, the spread in moments is very similar to the fixed-i^ ensemble. 

We now discuss the behavior of eigenstates in Hamiltonians with mixed K. In the locked- 
spin models, the K mixing can be induced by a tipped spin with 6 nonzero, or by a fully 
asymmetric rotor with the additional term eq. (j^) in the Hamiltonian. The effect of the 
tipped spin may be seen in the energy level plot in Fig. 3. The plot show the energies 
of the M = 1 states with spin angle tipped by 5° with respect to the symmetry axis. 
Several avoided level crossings may be seen. For example, the uppermost crossing involves a 
transition between states that have quantum numbers (JK) = (20) and {JK) = (1, —1) in 
zero field. We can see that the effect of the K-mixing is limited to regions very close to the 
level crossings; the spectrum is hardly changed otherwise. This means that the magnetic 
moments of the levels are unaffected except for special values of the field. The consequence 
is that the dispersion in magnetic moments is very similar to that of the fixed K ensemble. 
This may be seen in Fig. 4, where we show the magnetization distribution for a range 
of tipping angles. There is very little difference between the different Hamiltonians. The 
average magnetization from the diagonalized Hamiltonians with mixed K is also very close 
to the mixed- ensemble shown in Fig. 2. 

Returning to the spread in magnetization distribution, it appears that the locked-spin 
Hamiltonian cannot explain superparamagnetic behavior. The one caveat in this statement 
is that we have not considered fully asymmetric Hamiltonians. The perturbation eq. (9) does 
not produce a large mixing of the wave functions, and seems unlikely to have a significant 
effect on the widths. 



7 



IV. RESULTS: INTERMEDIATE COUPLING 



We now generalize from the locked-spin Hamiltonians to consider the spin as an inde- 
pendent degree of freedom. We shall treat the intermediate coupling spin-rotor Hamlitonian 
assuming that the inertia tensor is axially symmetric. As we saw in the last section, the 
breaking of K symmetry has no qualitative effect on the magnetic response, but only changes 
somewhat the magnitude of the susceptibility coefficient. We study the ensemble with a 
maximum orbital angular momentum of Jmax — 30 and the spin value S = 1. Despite the 
modest size of these numbers, this space requires diagonalization of about 4000 matrices of 
dimensionality ranging up to about 100. The ensembles are generated with a temperature 
of T = lOO/Ji as before. 

The first task is to see whether the locked-spin results can be reproduced with a spin S as 
small as 5" = 1. For this purpose we calculate the moment distribution taking a large value 
for the anisotropy coupling, A — b/iB. Typical results are shown in Fig. 5. In plotting the 
distribution, the magnetic moments have been scaled so that the maximum is ±1 in either 
model. We see that the two curves are practically identical, giving us confidence in using 
the Hamiltonian with a small value for S. The moment distribution in a weak coupling 
situation, A = O.Ol/^o-B, is shown in Fig. 6. Here the spin behaves as though it were 
uncoupled, so one sees the three M-states as in the atomic defiection. Notice however that 
the occupation probabilities are unequal. The probabilities are shifted at the level crossings, 
which become avoided crossings due to the (small) anisotropy term. At intermediate values 
of the coupling, the distribution changes smoothly from the three-peak structure in Fig. 6 
to the typical locked-spin shape of Fig. 5. In all cases the distribution remains broad. 

The field dependence of the magnetization shows a more complex behavior with inter- 
mediate coupling than we found for the locked spin. There are several distinct regions 
of behavior, depending on the relative sizes of the magnetic energy with respect to the 
anisotropy energy, and the rotational frequency with respect to the Larmor frequency. In 
extremely low fields, with magnetic energies less than the anisotropy energy, the spin should 
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behave in a similar way to the locked-spin model. When the magnetic field energy becomes 
larger than the anisotropy energy, the spin decouples from the rotational motion and we 
should see the individual Sz states as in Fig. 6. These two regimes are illustrated in Fig. 
7, showing the magnetization distribution for x = 0.005 and 0.02 for a model with weak 
coupling, A = O.OlkT. The distribution in the lowest field, shown by the solid line, is broad 
with a peak at ;U = 0. This is similar to the locked-spin behavior, ~ \og{^/^o). At 

X = 0.02, where the magnetic field is twice the anisotropy energy, the distribution has al- 
ready dissolved into the 3-peak structure of the decoupled S = 1 states. The behavior of 
the average magnetization as a function of x is shown in Fig. 8. In the extreme low- field 
limit, fioB < A, the magnetization is close to the weak-field Langevin function, eq. (3). 
The spin decouples from the rotational motion when B exceeds A, as already mentioned. 
From Fig. 8 it may be seen that the average magnetization increases rather slowly between 
yUo-B ~ O.OlkT and O.OSkT. Above that point the magnetization rises to even exceed the 
Langevin value. We believe that this change of behavior is associated with level crossings of 
the type {J Sz) ^ ( J ± 1 5*2 =F 1), corresponding to the matching of rotational and Larmor 
frequencies. A change of behavior at this matching point has been suggested previously 

It is interesting to examine the susceptibility in the third region in more detail. Here the 
magnetization is still far below the saturation condition fiQB ~ kT, and it is roughly linear 
in the field. We make a linear fit to the magnetization with the function, 

/i = afiQX, (12) 

The dependence of the fit coefficient a on magnetic field is plotted in Fig. 9. We see that the 
behavior like the the locked-spin model only if the anisotropy energy is substantially above 
kT. Below that, the coefficient exceeds the Langevin value, a = 1/3, as we saw on Fig. 7. 



Our findings for the magnetic response are quite different from those of ref. [jTT|, where 
a similar Hamiltonian was considered. In that work, the magnetization was found to be 
independent of the temperature but to depend strongly on the anisotropy energy. However, 
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these authors did not treat the adiabatic ensemble, but rather a mixed ensemble in which 
the magnetic field was turned on suddenly. 

We conclude this section with an estimate of the magnitude of A for iron clusters. The 
cubic anisotropy constant K for bulk iron is measured as ^ 4 x 10^ ergs/ cm^ = 5 x 10~^ 

eV/atom. For a cluster of 20 iron atoms this is 2.5 x 10~^ eV, which scaled to room 
temperature implies A ~ O.Ol/cT. Thus the weak coupling limit we discussed above should 
apply to moderate sized iron clusters at room temperature. A typical rotational frequency 
is given by (J^)^/^// ~ (A;T//)^/^. This has a magnitude of about 10~^ eY /h for a cluster 
of 20 iron atoms, which is the order of magnitude of the Larmor frequency in a 1 T field. 
Thus, the model predicts that light mass iron clusters would behave with the spin nearly 
decoupled from the rotational motion and with a fairly weak average susceptibility. This 
is certainly contrary to experiment, which has never demonstrated atomic-like behavior in 
clusters larger than diatomic. 



V. CONCLUSION 

Our study of the magnetic response of rotors coupled to spin is guided by the question, 
to what extent can a thermal ensemble of rotors play the role of the thermal bath that 
one assumes in the Langevin model? One's intuition may be that the rotor has far too few 
degrees of freedom to serve as a heat bath. Nevertheless, the low-field magnetic susceptibility 
of the thermal bath is reproduced on average by the deformed rotor with a locked and tipped 
spin axis. It is also obtained for intermediate coupling of spin to the rotor, except in a small 
interval of field strengths. 

However, superparamagnetism implies not only that the average moment is given by 
the Langevin formula, but all clusters should behave identically. Here the adiabatic rotor 
ensemble behaves quite differently, always having a large variance in the magnetization. 

Let us ask now, what would be needed in the spin-rotor Hamiltonian to achieve su- 
perparamagnetism? Clearly there would have to be much stronger K mixing than can be 
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generated by the spin-rotor Hamiltonian. Recall in Fig. 3, the energy levels have differ- 
ent slopes with large stretches of constant moment separated by small intervals of the level 
crossings. 

In the current idiom, we may speak of this as a weakly perturbed regular Hamiltonian. 
To achieve superparamagnetism, the Hamiltonian would have to be chaotic. Quantum 
mechanical chaos is a regularity and evenness in the energy spectra associated with the 
complete mixing of wave functions. In such a spectrum, the slopes of the levels would be 
similar and governed by the average behavior of the levels in that region of energy. When 



the Hamiltonian is chaotic, the slopes have universal correlation properties This 
behavior would require additional terms in the Hamiltonian eq.(^. A minimum condition 
for chaotic behavior is that the off-diagonal matrix elements exceed the level spacing of the 
regular part of the Hamiltonian ^M. 



It is expected that the coupling of the spin to vibrations could provide the thermal bath 



that would produce superparamagnetism |T9|,|20[]. However, it has not been demonstrated 



that this mechanism actually would work in finite-size clusters. The first issue is that vibra- 
tional transitions would have to be at the same energy as the spin flip in the external field. 
This is to allow off-diagonal matrix elements to nearby levels. The maximum vibrational 
frequency has the order of magnitude of the Debye energy huD] the spacing of vibrational 
levels AE^ih ~ huJo/<iN, where is the number of atoms in the cluster. In a 1 T field, the 
spin flip energy is 10~^ eV. With a Debye energy in iron of 6 x 10~^ eV, the vibrational level 
spacing is smaller than the spin- flip energy for clusters larger than = 20. Spin relaxation 
in the bulk takes place mainly through Raman scattering of phonons. In this case one is 
interested in the difference of two vibrational energies, and the vibrational level spacing 
again sets the scale. The next condition is that the average off-diagonal matrix element of 
the vibrational coupling should exceed the level spacing. This will always be satisfied for 



sufficiently large clusters, because V and AE vary as V ^ y (AE) as the size of system 
varies. This follows from the requirement that the relaxation time is independent of the size 
of the system for large A^. Unfortunately, the coupling is not well enough known at present 
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to make a numerical estimate. 
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FIGURES 

FIG. 1. Comparison of magnetization distributions of spherical clusters and deformed clusters, 

for a range of inertia ratios Ixjlz- 

FIG. 2. Comparison of average magnetization of mixed-ii' ensemble with fixed-i^ ensemble and 
with the Langevin function, eq. (1). 

FIG. 3. Energy levels for locked-spin Hamiltonian as a function of magnetic field, for low M = 1 
states. The spin is tipped by 5° with respect to the symmetry axis. 

FIG. 4. Magnetization distribution of deformed clusters, h/h = 2, with three different spin 
angles 6' = 0°,45°, and 90°. 

FIG. 5. Magnetization distribution of clusters in the anisotropic model, showing the strong 
coupling limit. Here the spin (s = 1) is coupled to the easy axes, with A = b/iB and x = 2. The 
distribution is close to the result of locked-spin model (dashed line). The distribution is plotted as 
a function of the rescaled magnetization JI. 

FIG. 6. Magnetization distribution for the weakly coupled anisotropic Hamiltonian. The spin 
value is 5 = 1, giving three peaks. 

FIG. 7. Magnetization distribution for weakly coupled spin to rotor, A = O.OlfcT. The solid 
curve shows the distribution in a very weak field, x = 0.005, and the dashed line shows the 
distribution in a field x = 0, somewhat stronger than the anisotropy energy. 

FIG. 8. Average magnetization of the anisotropic model for A = O.OlkT, as a function of 
Langevin parameter x = fiB/kT. 

FIG. 9. Fit of the low-field adiabatic susceptibility coefficient a in eq. (12) as a function of the 
anisotropy energy. The temperature is fixed at T = 50. 
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